Ab initio study of a mechanically gated molecule: Prom weak to strong correlation 
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The electronic spectrum of a chemically contacted molecule in the junction of a scanning tunneling 
microscope can be modified by tip retraction. We analyze this effect by a combination of density 
functional, many-body perturbation and numerical renormalization group theory, taking into ac- 
count both the non-locality and the dynamics of electronic correlation. Our findings, in particular 
the evolution from a broad quasiparticle resonance below to a narrow Kondo resonance at the Fermi 
energy, correspond to the experimental observations. 



Both the control of the geometric structure of a molec- 
ular junction and the systematic manipulation of its elec- 
tronic structure by an external parameter are of central 
importance for molecular electronics [![. So far, the best 
control over the geometric junction structure is reached 
in experiments based on scanning tunnelling microscopy 
(STM) Q, Q, because they allow the selection of an in- 
dividual molecule in a specific environment and its con- 
tacting with the STM tip at a defined position within the 
molecule (e.g. ref |3|). On the other hand, the most com- 
mon approach to tuning the electronic structure, electri- 
cal gating, is difficult to combine with STM. However, in 
a recent experimental study Q we contacted a surface- 
adsorbed molecule with an STM tip and peeled it off 
the surface by tip retraction, as shown schematically in 
Fig. [T] Spectroscopic data recorded during tip retrac- 
tion revealed a mechanical gating effect, in the sense that 
one of the molecular orbitals responds to the structural 
change and shifts with respect to the Fermi level (Ep) of 
the substrate, before becoming pinned at Ep (Fig. [2jl). 

In this letter we use this experiment, which was car- 
ried out on 3,4,9, 10-perylene-tetracarboxylic-dianhydride 
(PTCDA) adsorbed on the Ag(lll) surface, as a moti- 
vation for a theoretical study of the interplay between 
the geometric structure and the electronic spectrum of 
a molecular junction. The mechanical gating is subject 
to subtle details of electronic correlation, e.g. screen- 
ing of the intramolecular Coulomb repulsion by the elec- 
trodes, which is not described by conventional density- 
functional theory (DFT). We therefore propose the fol- 
lowing strategy: DFT addresses the atomistic details of 
the junction structure but does not provide reliable, well- 
founded spectral data. To evaluate the electronic spec- 
trum, we combine DFT with many-body perturbation 
theory (MBPT) [6] for non-local correlation and with 
numerical renormalization-group (NRG) theory Q for 
correlation dynamics beyond the mean-field level. We 
stress that only by combining DFT, MBPT and NRG a 
well-founded picture emerges, which can by systemati- 
cally compared to the data presented in Ref. [5| and the 
(more comprehensive) experiment in Fig. [2ji. 

We first discuss the geometric structure of the 
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FIG. 1: DFT-LDA data of tip/PTCDA/Ag(lll) junction, 
(a, b) Two representative configurations at tip-surface dis- 
tances of 2 = 7 A and z= 10 A, respectively. For better visi- 
bility only a section of the tip and the substrate are shown, 
(c) Bond energies of the PTCDA-tip and PTCDA-substrate 
bonds as a function of tip height z. At each z the molecule 
has been fully relaxed. The internal deformation energy of 
the PTCDA reflects the bending and twisting of its perylene 



tip/PTCDA/Ag(lll) junction in the framework of the 
local density approximation (LDA) to DFT Q. PTCDA 
molecules adsorb on the Ag(lll) surface in a flat-lying 
configuration Several DFT studies of the adsorp- 

tion have been published (for details, see Ref. [3| and 
references therein). Following experiment we place 
the tip, which in our simulation consists of ten Ag atoms 
in a pyramidal shape (IBj . above one of the carboxylic 
oxygen atoms and approach it to the surface. At a tip- 
surface distance of z = 6.2 A (6.7 A in experiment [5|) 
the oxygen atom jumps up and forms a covalent bond 
with the tip apex atom. The tip can then be moved up 
and down reversibly, forcing the oxygen atom and the 
attached section of the PTCDA molecule to follow, re- 
sulting in a peeling-like motion (Fig. QJi-b). Throughout 
the paper, z specifies the vertical distance between the 
tip apex atom and the uppermost surface layer. 
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Some representative DFT-LDA data of the structure 
and energetics of the junction are displayed in Fig. [1] 
Fig. [1] a and b shows two typical configurations. At a 
tip-surface distance of z = 7 A only a small section of the 
molecule is detached from the surface, while at z = 10 A 
about half of the molecule has lost contact with the sub- 
strate. The continuous nature of the peeling process as 
a function of z is reflected in the smooth behavior of the 
molecule-tip and molecule-surface bond energies (Fig.QJi) 
and in a smooth z-dependence of all intramolecular struc- 
tural parameters. 

The most interesting feature of the electronic structure 
of the junction is the PTCDA LUMO (lowest unoccupied 
molecular orbital), because this is the orbital for which 
the gating effect as a function of the external parameter 
z is observed In the gas phase, the LUMO is found 
-2eV above the Ag(lll) Fermi level (Fig. HJi). Upon 
adsorption, the LUMO is lowered in energy below the 
Ag(lll) Fermi level due to the metallic polarizability of 
the substrate 111- 21 1 and thus becomes partially occu- 
pied with about 1.8 electrons (counterbalanced by back- 
donation from other molecular states to the surface). 

The partial occupation of the LUMO (which changes 
as a function of z; see below) causes the high sensitiv- 
ity of the LUMO spectrum to the junction structure: In 
experiment, for z below 8.2 A the LUMO approaches 
the Fermi energy at a steep slope of ~0.2-0.3 cV/A and 
turns into a sharp resonance as soon as it reaches Ep at 
z=8.2 A. At z =9.7 A, the FWHM of the LUMO peak 
is 14 meV. We note that although a mechanical gating 
effect is observed already in DFT-LDA (top curves in 
Fig. [2j>c and open circles in Fig. [2ji), there is a large 
discrepancy with the experimental data: The DFT-LDA 
LUMO reaches the Fermi level far too late (at 12 A), and 
it does not sharpen as much as in experiment. Because of 
this sharpening, the experimental resonance at the Fermi 
level has been discussed in terms of dynamic correlations 
(Kondo effect) @. 

A systematic approach to the Kondo effect is achieved 
by taking the single-particle mean-field spectrum of the 
LUMO (at each z) as a starting point for a NRG sim- 
ulation (see below for details). However, the DFT-LDA 
spectrum should not be used for this purpose since it 
lacks physical significance, being essentially an auxiliary 
quantity within a ground-state total-energy approach. A 
proper single-particle mean-field spectrum to be used as 
input to NRG must refer to electronic excitations (i.e., 
removal or addition of one electron) Q. In particular, 
these excitations are subject to non-local correlation ef- 
fects that are not included in DFT-LDA. 

Mean-field spectra which take non-local correlation 
into account can be addressed within the standard GW 
approximation of MBPT [§}. In the present case, how- 
ever, a full GW calculation is too demanding. We 
therefore employ a simplified, perturbative LDA+GtflU 
approach which yields reliable QP energies E Gdw := 
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FIG. 2: (a) HOMO and LUMO energies of gas-phase PTCDA 
and PTCDA/ Ag(lll) (peak of the PDOSpliwithout STM 
tip). Experimental data are from Ref. \YA, [la ] (gas phase) 
and Ref. (on Ag(lll)). The two energy scales differ by 
Ef — -Evac — Ew (Ew — work function of Ag(lll)). (b,c) 
LUMO spectra for the tip/PTCDA/Ag(lll) junction [r| at 
tip heights of z = 8 A and z = 11 A. The experimental spectra 
are cuts through the dl/dV map of panel d, at z-values of 7.5 
A and 9 A. (d) LUMO peak positions and experimental dl/dV 
color map, ranging from 0Go (light blue) to 0.10 Go (yellow). 
Experiments were performed by approach-retraction cycles of 
the tip, with bias voltages increasing in steps of 1 mV. dl /dV 
was detected with a lock-in amplifier, modulation amplitude 
4mV and modulation frequency 2.9 kHz. For bias voltages 
exceeding ±100 meV the tip-PTCDA contact becomes unsta- 
ble. The experimental data point at -0.45 eV is an estimate 
from Ref. @ (-0.35 eV with the tip attached to the PTCDA 
monolayer) and Ref. [12 (difference -0.1 eV between single- 
molecule and monolayer spectra). 
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A9 dw on top of DFT-LDA from a QP Hamilto- 
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= H 



LDA 



iG(W - W metal ) 



(1) 



in which the QP self-energy correction A Gdw := iG(W — 
Wmetal) results from the difference between the correctly 
screened Coulomb interaction (W) and an interaction 



Wmetal from hypothetical metallic screening (see Ref. [21 1 
for details). The resulting QP corrections to DFT-LDA 
differ from standard GW QP corrections by less than 20 
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FIG. 3: Properties of the PTCDA LUMO at varying tip- 
surface distances z. (a) Mean-field QP correction (from 
LDA+GdW) of the LUMO for various occupancies, i.e. (from 
top to bottom) for the state being empty, being partially oc- 
cupied (at n\_DA+Gdw), an d being filled, (b) Intra-state in- 
teraction U (from Eq. ©). (c) Occupancy of the LUMO, 
resulting from LDA and LDA+GdW . Note that the number 
of electrons in the state is 2-n. (d) Resulting Kondo temper- 
ature. 



% (2lH2J|. For gas-phase PTCDA, LDA+GdW can eas- 



ily be checked against a standard GW calculation and 
against the experiment. Fig. [5^ shows that the GdW 
shifts (A Gdvv ' = -1.8cV for the HOMO and +1.7 eV for 
the LUMO) agree well with standard GW results (both 
with our own and with those of Ref. [17J), in particular 
for the LUMO, and with experiment [17|,|l8|, establishing 
LDA+GdW as a suitable method. 

When the molecule is adsorbed on the surface, three 
features show up in A Gdw and E Gdw : (i) Due to the 
metallic nature of the substrate, which screens the in- 
teraction inside the molecule non- locally HHH, the 



GdW corrections are substantially reduced as compared 
to gas-phase PTCDA (25|. As shown in Fig. |3k, the QP 
corrections are below 0.5 eV for z = 5 A instead of al- 
most 2eV for the free molecule. Note, however, that as 
the molecule is removed from the metallic substrate, its 
screening becomes less metallic, and both W and A Gdw 
grow again with increasing z. (ii) The GdW correction 
of the LUMO changes sign as a function of its occupa- 
tion. If the state were empty (occupied) it would observe 
a GdW correction of A G £Jf >0 (A G if <0), as shown in 
Fig. [3^,. At partial occupation 0<n<l (through charge 
transfer from the metal) we obtain an n-dependent QP 
shift A Gdw {n) = (1 - n)A G i™ + nA Gd ™ . (hi) Since 
the QP correction A Gdw (n) shifts the whole LUMO 
spectrum f(E) rigidly with respect to the Fermi level, 
thereby changing its occupation by An, the classical 
(screened) Coulomb energy changes by 2(An)U (the fac- 



tor of 2 results from spin degeneracy), which must be 
added to A Gdw , yielding at effective GdW correction of 
A Gdw {n) = (l~n)A Gd ^+nA G £f+2(n~n hDA )U. Since 
the QP shift in turn determines n = /f£ f Gdw {E)dE via 
f Gdw {E) = f LBA (E - A Gdw {n)), n and A Gdw must be 
determined self-consistently for each z (Fig. [3£i,c). At 
small z we obtain a negative QP shift and a slight in- 
crease of n (as compared to jt-lda) while at large z the 
QP shift approaches zero. A key ingredient into our cal- 
culations is the intra-state interaction energy U . It is 
given (for each z) by 



U 



\i> W Mo(r)\ 2 W(r,r')\^ hVMO (r')\ 2 d 3 rd 3 r' . (2) 



Fig. reveals that as a consequence of efficient metallic 
screening the U parameter is rather small (below 1 eV) 
for z<8A, but increases for larger z. For gas-phase 
PTCDA it amounts to 3.0 eV. 

Results for the effective QP shift A G g w , the U param- 
eter, and the self-consistent occupation n are displayed 
in Fig. [3^,-c. Although the QP correction turns out to 
be small, there is a clear trend from A G ^ W = — 0.17eV 
for z = 5 A to A Gdw ~ eV for z = 12 A. This leads to an 
increased slope of the LDA+GdW QP peak position for 
z between 7 and 10 A (Fig. Hi). 

We finally extend our LDA+GdW calculation by in- 
cluding dynamical correlation in the NRG approach. Us- 
ing the calculated values of U(z), Eq. ([2]), we are able 
to extract the bare LUMO level positions eo{z) and cou- 
pling functions T(E, z) from the LDA+GdW spectra at 
each z, by equating 



f UdW (E,z) 



1 



7T E-e (z)-nU(z)-T(E,z) 



(3) 



Together with U(z), the so obtained eo(z), T(E, z) enter 
the NRG calculation 0, [26[ which yields the many-body 
LUMO spectra displayed in Fig.[2j>c (27|. 

The z-dependent renormalizcd peak position of the 
LDA+GdW+NKG data (full squares) can now be com- 
pared to the experimental data as shown in Fig. [2ji. As 
in experiment, the calculation shows a fast-rising QP 
peak (slope 0.15 eV/A for z between 8 and 10 A), which 
also becomes pinned at Ep as soon as it reaches the 
Fermi energy (at z w 10 A). In both experiment and 
LDA+GgW+NRG, the pinned resonance is substantially 
sharper than the QP peak that moves up to Ep (see Fig. 
[2]c). Since this sharpening to 14 mcV (experiment at 
9.7A) or 15 meV (theory at 11 A) is neither observed in 
LDA nor in LDA+GdW, we can ascribe it to a dynamical 
correlation effect. 

Note that for z<&A the LDA+GdW and 
LDA+GdW +NRG spectra are very similar (Fig. HJd). 
This is to be expected, since a nearly fully occupied state 
below Ep shows very weak dynamical correlation. In 
contrast, for z>9A the difference between LDA+GdW 
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and LDA+GdW+NRG becomes dramatic. While 
LDk+GdW still yields a single LUMO peak (Fig. [2p), 
the LDA+GilU+NRG spectrum shows a three-peak 
structure, with two single-particle side bands (at -0.9 
eV and +0.3 cV for 2=11 A, separated by ~J7) and a 
Kondo resonance at Ep (Fig.^). The Kondo resonance 
is also observed in experiment (bottom of Fig. [2J; and 
sharp horizontal line in Fig. [2jl). 

In Fig. |3ji the evolution of the Kondo temperature 
Tx(z) is plotted. Tk has been determined from the NRG 
calculation as the temperature at which 40 % of the local 
moment in the LUMO is screened (Wilson criterion). For 
z<9A T K is a few thousand K, because in this z-range 
U and r are of similar size, and hence there is no clear 
separation between spin and charge fluctuation energy 
scales. In other words, in this regime Tk is equivalent 
to the charge fluctuation scale T. Physically, this means 
that screening of the small residual moment that is left 
by the total occupancy of 1.8 electrons in the LUMO is 
accomplished by uncorrelatcd charge fluctuations in and 
out of the level. At z = 12 A T K has dropped to T K < 1 K, 
caused by the increase of the ratio U/T which leads to 
the clear separation of spin and charge energy scales. Ac- 
cordingly, Tk is now a Kondo temperature in the proper 
sense, related to the energy scale of correlated screening 
of a fixed moment in the LUMO by virtual spin fluc- 
tuations. In experiment, the Kondo resonance can only 
be followed up to z w 10 A because of decreasing signal 
strength, but at 9.7 A the Kondo resonance has a FWHM 
of 14 meV, in good agreement with a predicted Kondo 
temperature of the order 100 K (Fig. |3ji). 

One remaining difference between our 
LDA+GdVF+NRG data and experiment is the ap- 
position at which the QP peak reaches Ep and turns 
into the resonance. This could be due to the under- 
estimation of the Ag(lll) work function within the 
underlying DFT-LDA (^4.6 eV instead of the experi- 
mental value of 4.8 eV), which persists throughout our 
calculations. Taking this into account, the Fermi level in 
all theoretical spectra should be lowered by 0.2 eV, i.e. 
all molecular QP features in Fig. [2] would be at higher 
energies by 0.2 eV. The LDA+GdIT+NRG peak would 
then reach the Fermi level at about z ~ 8 A. 

In conclusion, we have presented an approach 
to calculating electronic spectra that systematically 
includes non-local and dynamic correlation effects. 
It correctly predicts the mechanical gating of the 
tip/PTCDA/Ag(lll) junction. Because of its compu- 
tational efficiency, this state-of-the-art approach to elec- 
tronic spectra should be widely applicable. 
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